NMR relaxation rate and dynamical structure factors in nematic and multipolar 
liquids of frustrated spin chains under magnetic fields 
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Recently, it has been shown that spin nematic (quadrupolar) or higher multipolar correlation func- 
tions exhibit a quasi long-range order in the wide region of the field-induced Tomonaga-Luttinger- 
liquid (TLL) phase in spin-i zigzag chains. In this Rapid Communication, we point out that the 
temperature dependence of the NMR relaxation rate l/Ti in these multipolar TLLs is qualitatively 
different from that in more conventional TLLs of one-dimensional quantum magnets (e.g., the spin-i 
Heisenberg chain); l/Ti decreases with lowering temperature in multipolar TLL. We also discuss 
low-energy features in spin dynamical structure factors which are characteristic of the multipolar 
TLL phases. 

PACS numbers: 76.60.-k,75.40.Gb,75.10.Jm,75.10.Pq 



Magnetic states with an order parameter defined by 
a product of multiple spins, such as nematic, vector 
chiral, and scalar chiral orders, have attracted much 
attention. Spin nematic (quadrupolar) ordered phases 
have been recently shown to appear in frustrated fer- 
romagnets, such as ferromagnets with competing anti- 
ferromagnetic (AF) interactionsii^i^ and magnets with 
multi-spin-exchange couplings.— >^>^ A triatic (octupolar) 
ordered phase was also found in the triangular lattice 
multiple-spin exchange model with ferromagnetic (FM) 
dominant coupling)^ These nematic and triatic ordered 
states can be regarded as Bose condensed states of bound 
two-magnonsii^ and bound three-magnons;^ respectively, 
and their order parameters are given by Sj'S^ or S~ Sj^ , 
and S^S'^S^ or Sj^ Sf' . 

Recent extensive studies^'^'^'^iiS have shown that a 
series of similar multipolar phases appear in the one- 
dimensional (ID) spin-i Heisenberg model with FM 
nearest-neighbor exchange Ji and competing AF next- 
nearest-neighbor exchange J2 in applied field H, whose 
Hamiltonian is 



(1) 



Here Sj is the spin-i operator on jth site, Ji < 
and J2 > 0, and H is the external magnetic field 
in the z direction. This simple frustrated spin chain 
is a minimal model of frustrated ferromagnets and is 
thought to describe magnetism in quasi-lD edge-sharing 
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such as Rb2Cu2Mo30i2, NaCu202, 



LiCuV04, and LiCu202. 

Hikihara et al^ and Sudan et ali^ showed that 
the ground state of Hamiltonian ([T]) has field-induced 
Tomonaga-Luttinger (TL) liquid phases in which spin 
multipolar correlations are quasi-long-range ordered 
while the transverse spin correlation is short-ranged. 
This result can be easily understood in the large- 
magnetization regime, where p magnons form a bound 
state (p > 1).— 1^ A gas of bound p magnons acquires 
off-diagonal quasi-long-range order, with the order pa- 



rameter being the effective hard-core boson creation op- 
erator ri/=i •■■ p^j+i- original spin language this 
off-diagonal correlation is the multipolar spin correlation 
characterizing a nematic (p = 2), an octupolar (p = 3), or 
a hexadecapolar {p = 4) phase. Numerical studies^i^iii^ 
found p = 2 for -2.7 < J1/J2 < 0, p = 3 for -3.5 < 
J1/J2 < -2.7, and p = 4 for -3.76 < J1/J2 < -3.5, 
at the saturation field. At lower magnetic fields these 
phases cross over to spin-density wave SDWp phases in 
which the density correlation of bound p magnons, i.e., 
the longitudinal SDW, becomes stronger than the multi- 
polar correlation. Incidentally, a SDW2 phase is present 
also in the case of AF Ji > 

However, it will be difficult to obtain direct experi- 
mental evidence for the multipolar spin orders, as it re- 
quires probing four- or more-spin correlation functions 
with high accuracy. Standard experimental probes, such 
as neutron scattering or magnetic resonance, measure 
only two-spin correlations. Furthermore, the multipo- 
lar TL liquids have a gapless spectrum and a smooth 
magnetization curve. Thus, if one only measures their 
static, thermodynamic quantities (uniform susceptibility, 
specific heat, entropy, etc.), it is hard to distinguish the 
multipolar TL liquids from conventional TL liquids. Ex- 
perimental schemes for identifying multipolar spin orders 
are therefore called for. 

In this Rapid Communcation, we propose that NMR 
measurements can capture signatures (albeit indirect) of 
the multipolar TL liquids. We show that, in the TL 
liquids with a dominant multipolar spin correlation, the 
NMR relaxation rate 1/Ti decreases as temperature T 
is lowered [sec Fig. [IJa)]. This temperature dependence 
of 1/Ti is opposite to that in conventional TL liquids 
(and in SDWp phases), where it is always diverging as 
T — > 0, in magnetic fieldiii We also point out that the 
spin dynamical structure factors exhibit features which 
are very characteristic of the multipolar TL liquids. 

Let us begin with a brief review of the effective theo- 
riea^i^ for the multipolar TL liquids, which allow us to 
find low-energy behavior of spin correlation functions ob- 
served in NMR and in dynamical structure factors. In 
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FIG. 1: (color online) Temperature and field dependence of 
the NMR relaxation rate l/Ti in (a) multipolar phases in 
spin-i frustrated zigzag chains and in (b) standard TL liq- 
uids (e.g., the AF Heisenberg chain in magnetic field). The 
parameter K controls correlations in TL liquids (see the text). 



the weak Ji limit, the Abelian bosonization methodic is 
useful. It leads to nematic {p = 2) and SDW2 phases^i^ 
and a vector chiral ordered phaseJ^ The effective Hamil- 
tonian for the nematic {p — ■^"'^ QriAAr_ ^koooofLI^ 
written as 



2) and SDW2 phasesS^iS is 



+ 5sin(7rM)sin(\/87r(?!)_ +7rA/)|, (2) 



where x = 2j (the lattice spacing is set equal to unity), 
i4>±i(^±) is a pair of dual scalar fields satisfying the com- 
mutator [(j)^{x),dy9^{y)] = i6^^J{x - y), M = {S^, 
5 (X Ji, and K± and v± arc, respectively, the TL-liquid 
parameter and the velocity of the {(f>±,0±) sector. The 
{4>+, 9+) sector is a gaplcss TL liquid, while the (0_, 6^) 
sector has a gapful spectrum because the field c/)- is 
pinned at a value minimizing the potential energy of the 
g term. Using Hamiltonian ([2]) , one can evaluate the low- 
energy and long-distance behaviors of several correlation 
functions. The imaginary-time (r) spin and nematic cor- 
relations at zero temperature T = are calculated as 



(5|(r)5o^(0))=A/ 
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(5;(r)5o-(0)) = C2COS (f) , ,vW)U 11/2 + 



(3a) 



(3b) 



(5+(r)5++i(r)5o-(0)5r(0)) = ^3^^^^ + • • • , (3c) 

where z± — j ~~ iv±T, and C„ are nonuniversal positive 
constants. The exponential decay of the transverse spin 
correlation in Eq. (|3b[) is qualitatively different from a 
power-law decay form in ordinary TL liquids (e.g., the 
spin-i AF chain and ladder in magnetic field). The cor- 
relation length ^ is inversely proportional to the gap of 
the {(/)-, 6 _) sector. 

More generally, the TL-liquid behavior in all the multi- 
polar and SDWp phases {p > 2) can be understood from 
a hard-core Bose gas picture of bound p magnons, when 
the nearest-neighbor coupling Ji is fcrromagnetici^ Be- 
low the saturation field a (dilute) Bose gas of bound 
p magnons form a TL liquid with off-diagonal quasi- 
long-range order, i.e., pth multipolar TL liquid. In 
this picture one may replace the pth multipolar opera- 
tor S~^j^S~^_2 ■ ' ' ^j+p magnon density ■7; ^ Sj with a 



creation operator of a hard-core boson (— l)^6j and bo- 
son density pb^jbj, respectively. Here the staggered fac- 
tor (—1)-' represents the total momentum A: = tt of the 
lowest-energy bound states. The hydrodynamic theory 
for the bosonic TL liquid has the same form as the free 
boson Hamiltonian of the {(f>+, 9+) sector in Eq. ([2|). The 
effective theory gives the following longitudinal spin and 
the multipolar correlation functions at T = 0: 



p^K f 1 
4^ 
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n=l 



(1 - 2M) 

>|1/(2A') 



(4a) 
(4b) 



where z = j — ivr, K is the TL-liquid parameter for the 
hard-core bosons. While we cannot evaluate the trans- 
verse spin correlations within this boson picture, they 
must decay exponentially cx cxp(— |z|/^) as it is neces- 
sary to break a magnon bound state in order to create 
an excitation with AS*^ = ±1. In the nematic case of 
p = 2 and Ji < 0, Eq. (|4]) coincides with Eqs. ((3a|) and 



if we set = 2K and 



V. Near the saturation 



field where the density of magnons vanishes, the value of 
K approaches unity, i.e., that of the ID free fermions. 



3 



Indeed, the numerical calculations in Ref. have shown 
that K monotonically increases from about 1 /4 to unity 
with the increase in the magnetization M . This means 
that the multipolar correlation (|4b[) is strongest in the 
high-field regime {2K > 1), whereas the SDW correla- 
tion ()4ap becomes most dominant in the low-field regime 
{2K < 1). This property is important in the following 
discussion on the NMR relaxation rate. We note that at 
p = I, Eq. ([4]) reproduces the spin correlations in the TL- 
liquid phase of, e.g., the spin-i AF chains under magnetic 
fields 

The temperature dependence of the NMR relaxation 
rate 1/Ti in the multipolar TL liquids can be derived 
from the above asymptotic forms of correlation functions. 
The perturbation theory in hyperfine interaction between 
nuclear and electron spins obtains 1 /Ti a o^"^'^" 



1 



E 



[S+-{k,uj) + S-+{k,uj)] 



\Al\ S'^'ik^iu) 



(5) 



where to is the nuclear resonance frequency, A'jl 
are the hyperfine form factors, and S°'^{k,Lu) ~ 
Ej-e"'''^ /^d<e*'^*(S'j'(<)S'^(0)) is the spin dynamical 
structure factor [t — —it is the real time) at tempera- 
ture T. Since lu is generally much smaller than the en- 
ergy scale of spin exchange interactions, we may take 
the limit uj/T +0. Moreover, the k dependence of 
A'j^ is usually weak due to the locality of the nucleus- 
electron interaction. Hence, the T dependence of 1/Ti 
can be obtained by evaluating the local susceptibility. 

The local susceptibility at finite temperatures can be 
readily obtained from the correlation functions ^ and 
(HI) through the standard procedureJ^ii^ Substituting 
them into Eq. (O, we obtain 1/Ti for the multipolar TL 
liquids in the form 



1/Ti 



DlT 



(6) 



The two leading terms come from the second and third 
terms, respectively, of the longitudinal spin correlation, 
Eqs. ([5a|) or Pa|) . The coefficients d\ and d\ are inde- 
pendent of temperature in the regime uj <^ T l'^i.2|- 
In Eq. ^ we have omitted contributions from the trans- 
verse spin correlations which are exponentially small, 
g-A/T jg proportional to the spin gap), at 

low temperatures T ^ \ Ji.2\- When K < 1, the second 
term in Eq. ^ gives the leading contribution in the low- 
temperature limit. Similarly, the known T dependence 
of 1/Ti in spin-i AF chains under a magnetic field is 
obtained from Eq. ^ with p = 1 in the fori 



,17.18 



1/Ti = E^T- 



\rp2K-l 



^±j.l/{2K)-l 



(7) 



where the terms oc eI and E^ are derived from the longi- 
tudinal and the transverse spin correlations, respectively. 



Equation ^ commonly holds in TL-liquid phases of ID 
magnets such as AF spin chains and ladders in magnetic 
field.ii 

Comparison of Eqs. ([6|) and ^ tells us an important 
feature of the NMR relaxation rate in the multipolar TL 
liquids. As we noted above, the parameter K in the mul- 
tipolar phases of Hamiltonian ^ with the FM coupling 
Ji < 0, is an increasing function of H and approaches 
unity at the saturation fieldi^ The monotonic magnetic- 
field dependence of K presumably holds for other mul- 
tipolar TL liquids as well, at least for spin-i AF spin 
systems. Equation ^ then implies that 1/Ti decreases 
with lowering temperature in the high-field multipolar 
phase {2K > 1)^ while it shows diverging behavior in 
the low-field SDWp region {2K < 1); see Fig. [H This 
behavior is totally different from that of conventional TL 
liquids like AF spin chains under magnetic field [Eq. ([7])], 
in which 1/Ti always diverges in the low-temperature 
limit, irrespective of the value of K (the case of K = 1/2 
is specia l^^'^'^ ). We emphasize that this difference in the 
T dependence of 1/Ti between multipolar and conven- 
tional TL liquids, shown in Figs. [T] (a) and (b), can be 
taken as a pronounced signature of ID spin-^ multipolar 
TL-liquid phases. The decay of 1/Ti with lowering tem- 
perature in the multipolar liquid phases is due to both 
the absence of gapless modes in {k, us) and the weak 
singularity at w = in S'^^{k,uj) [see Eqs. (|8]) and 
We also note that NMR experiments cannot distinguish 
a SDWp region from ordinary TL liquids because they 
both show divergent behavior of 1/Ti as T — > 0. 

Next we discuss the spin dynamical structure factors 
S°'^{k,u!) at T = in the multipolar phases. The sup- 
port of S°'^{k,u!) tells us which excitations in the {k,Lu) 
space contribute to inelastic neutron scattering. The 
low-energy parts of ^"^(/c,^) are obtained from Fourier 
transform^! of correlation functions ^ and ^ . For the 
nematic and SDW2 phases in FM Ji < 0, wc find 

S''{k^0,uj)^4K\k\6{uj-v\k\), (8a) 

s-(k^±k.,.)= ;f;t::^':iu^ ^ (8b) 



[uj'^ +v'^{kTk2) 



^ ^'^^"/''")^ [c.-e(fcTV2)]^'V(4K,) -(8^) 

where fcs = 7r(l - 2M)/2, t{k) = {v^_P + A^^^, 6s(w) 
is a unit step function, and c^'^ are positive numerical 



constants. The (5-function peak in Eq. ((8a|) will have 
a finite width when the nonlinearity of the low-energy 
dispersion is included. In the SDW2 phase with AF Ji > 
0, in Eq. jSb)) should be replaced with k2 = 7r(l -|- 
2M)/2. The longitudinal part S''^{k,uj) in the higher 
multipolar and SDWp phases {p > 2) is also obtained 
from Eq. (Pa|) as 



S'^'ik - 0, w) = p'^K\k\5{uj - v\k\), (9a) 
c^6s(w — v\k =F kJ) 
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where 7 = + 



and c^j*- is a positive constant. 



These resuhs are depicted in Fig. O The gapless exci- 
tations giving dominant contribution to S^''^{k,uj) in the 
pth multipolar or the SDWp TL liquid are located at 
k = ±(1 — 2M)t:/p, when Ji < 0. These wave numbers, 
inversely proportional to the number p of magnons form- 
ing a bound state, are equal to the "2kp^^ of the hard-core 
Bose liquid of bound p magnons^ (note that fermions and 
hard-core bosons are equivalent in one dimension). The 
result for the ordinary TL liquid (e.g., the AF Heisenberg 
chain) corresponds to the case p = 1 , or the limit Ji — > 
(the lattice unit equals two in this case). Furthermore, 
one can discriminate between the SDW2 phases in Ji < 
and in Ji > by observing the shift of the gapless points 
from k = 7r/2 in S^^{k,uj). Another manifest difference 
between multipolar and ordinary TL liquids is that the 
transverse component S'"' {k,uj) has a gap in the mul- 
tipolar phases, while that of the ordinary TL liquids is 
gapless. These features in S°'^{k,uj) can be employed 
as definite signatures of the multipolar and the SDWp 
phases in the model ([T]). 



FIG. 2: (color online) Low-energy relevant parts of S^^{k,uj) 
and (fc, ui) in the multipolar and SDWp phases of the 
spin-i frustrated zigzag chain in magnetic field [(a), (b), (c), 
and (e)] and in the TL liquid of spin-i AF chains in magnetic 
field [(d) and (f)]. 



where kp = 7r(l — 2M)/p. For comparison. S°'^{k,uj) 
in the standard TL liquids, e.g., spin-i AF Heisenberg 
chains under a magnetic field, have the formic 



5^^(fc-0,t^) =K\k\5{Lu~v\k\), 

cf 9s(w — v\k =F fell 



5^^(fc-±A:i,tj) 
5+-(fc~±7r,w) 
S+~{k - ±2t:M,lu) 



(10a) 
(10b) 



[uj-i +v^{kTki)^]^-^'' 

[t^2 +^,2(fc^7r)2]i-i/(4A-) 

dj^Osi^ - v\kT2nM\) 

[Lu^v{kTlTrM)]^ 
X [u^ -v^ikTlT^Mfy , (lOd) 



To conclude, we have studied dynamical response of 
the multipolar TL liquids in the spin-i frustrated zigzag 
chains in applied magnetic field. The NMR relaxation 
rate 1/Ti in the multipolar TL liquids shows algebraic 
decay with lowering temperature, which is distinct from 
the diverging behavior in conventional TL liquids like the 
spin-i AF chains (see Fig. [1]). Furthermore the wave- 
number and the magnetization dependence of the gap- 
less modes in the dynamical structure factors S"^{k,uj) 
can provide us with clear evidence for the presence of 
the multipolar liquids as well as SDWp regions of bound 
magnons. Our arguments are also applicable to multi- 
polar phases in higher-spin chains, where K can become 
larger than imity due to soft-core repulsion of bosons. 
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